Interplay of couplings between antiferrodistortive, ferroelectric, and strain degrees of 
freedom in monodomain PbTiOs/SrTiOa superlattices 
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We report first-principles calculations on the coupling between epitaxial strain, polarization, and 
oxygen octahedra rotations in monodomain (PbTi03)n/(SrTi03)n superlattices. We show how the 
interplay between (i) the epitaxial strain and (ii) the electrostatic conditions, can be used to control 
the orientation of the main axis of the system. The electrostatic constrains at the interface facilitate 
the rotation of the polarization and, as a consequence, we predict large piezoelectric responses 
at epitaxial strains smaller than those that would be required considering only strain effects. In 
addition, ferroelectric (FE) and antiferrodistortive (AFD) modes are strongly coupled. Usual steric 
arguments cannot explain this coupling and a covalent model is proposed to account for it. The 
energy gain due to the FE-AFD coupling decreases with the periodicity of the superlattice, becoming 
negligible for n > 3. 
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Superlattices composed of thin layers of ferroelectric 
and paraelectric (or incipient ferroelectric) ABO3 per- 
ovskites oxides have attracted a lot of interest during the 
last few years [iHl]. The fascination for these layered 
systems comes from the fact that the properties of epi- 
taxial structures, made by stacking different perovskites, 
are not a simple combination of the properties of the con- 
stituent materials, but exotic phenomena might emerge 
that fully rely on interfacial effects. 

ABO3 perovskite oxides present different phase transi- 
tion sequences and ground states (GS) involving, among 
others, zone-center ferroelectric (FE) distortions, charac- 
terized by the opposite motion of the cations with respect 
the O cage, and non-polar zone-boundary antiferrodis- 
tortive (AFD) modes, which consist on rotation and tilt- 
ing of the oxygen octahedra surrounding the B-cation Q . 
But polar and nonpolar instabilities often compete and 
tend to suppress each other, so one of these distortions 
dominates over the others and is the only one that ap- 
pears in the GS structure. However, the balance is ex- 
tremely delicate and can be tuned by external electrical 
and strain [gl fields, or by changing the chemical envi- 
ronment through the use of different materials and peri- 
odicities in the stack Q- 

One of the most studied systems in the recent lit- 
erature, the (PbTi03)„/(SrTi03)„ (PTO/STO) super- 
lattice, constitutes a good example where the balance 
between different instabilities have been observed to be 
strongly tunable. It was theoretically predicted and ex- 
perimentally observed that the polarization, tetragonal- 
ity and phase transition temperature of the system can 
be monitored with the number of PTO layers, m, de- 
creasing monotonically when the PTO volume fraction is 
reduced [1,0. However, in the limit of ultrashort peri- 
ods, PTO/STO superlattices exhibit an unexpected re- 



covery of ferroelectricity that cannot be accounted for 
by simple electrostatic considerations alone [l^l- In this 
milestone work, the GS of the system was described as 
not purely ferroelectric, but involving a trilinear coupling 
term between two AFD modes, that correspond to in- 
phase (AFDzi) and out-of-phase (AFD20) rotations of 
the oxygen octahedra around the z axis, that induce a 
polar FE distortion (Vz) in a way compatible with hy- 
brid improper ferroelectricity. 

Despite all the previous efforts, two problems remain 
virtually unexplored. First, although many works have 
dealt independently with the FE-strain and AFD-strain 
couplings (see Ref. @ and [111, respectively) the influ- 
ence of the direction and magnitude of local polariza- 
tion on the rotation of the oxygen octahedra has not re- 
ceived the same attention. The control over the bond an- 
gles through the polarization orientation and magnitude 
could open new routes to generate or enhance magne- 
toelectric couplings [l2| . Second, first-principles simula- 
tions addressing the transition between the improper and 
normal ferroelectric regimes as a function of the periodic- 
ity of the superlattice are, up to our knowledge, missing 
in the literature. 

In this Letter we theoretically predict large mixed FE- 
AFD-strain couplings in monodomain PTO/STO super- 
lattices. As a result of these interplays, the phase dia- 
gram is much richer than originally assumed [loj . with 
rotation of the polarization away from the superlattice 
normal, and a strong coupling of the AFD modes with 
the magnitude and direction of the FE mode. We study 
the physical origin of this phase diagram and the energet- 
ics of the FE and AFD contributions and their coupling 
with increasing periodicity. 

For this study we perform first principles simulations 
of (PTO)„/(STO)„ superlattices, within the local den- 
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FIG. 1. (color online) (a) Squematic representation of a 
(2/2) PTO/STO superlattice. TiOe octahedra are labeled 
according to the chemical identity of the first two neighbor AO 
planes, and the direction of the polarization, (b) Definition of 
the angles of rotation around the z axis, (p, and tilting around 
an axis in the {x,y) plane, 0, of octahedra. 



sity approximation (LDA) using the Siesta code [13| . 
Very accurate computations are required since the rel- 
evant differences in energies are 7 orders of magnitude 
smaller than the absolute value of the energy. Real space 
integrations arc computed in a uniform grid, with an 
equivalent plane-wave cutoff of 1200 Ry. For the Bril- 
louin zone integrations we use a Monkhorst-Pack sam- 
pling equivalent tol2xl2xl2ina five atom perovskite 
unit cell. Details on the norm-conserving pseudopoten- 
tials and the basis set used can be found in Ref. 
The superlattices are simulated by means of a supercell 
approach, where we repeat periodically in space a basic 
unit cell, that is built stacking alternating n-unit-cells- 
thick layers of PTO and STO along the [001] direction 
for a global periodicity of {n/n) [Fig. [Ija)]. This struc- 
ture leads to naturally classify the TiOg octahedra into 
four different types (PTO, STO, P+ and P") depending 
on the top/bottom AO layer and the direction of polar- 
ization [see Fig. [TJa)]. In-plane lattice vectors are dou- 
bled to account for the condensation of AFD instabilities. 
With the (2 x 2) in-plane periodicity, TiOg octahedra are 
allowed both to rotate an angle </> around the z-axis or 
to tilt an angle 9 around an axis contained in the (a;,y) 
plane [Fig. [ijb)]. The mechanical boundary conditions 
imposed by the substrate are implicitly treated by fix- 
ing the in-plane lattice constant, a\\. The use of periodic 
boundary conditions imposes short-circuit electrical con- 
dition across the whole unit cell. 

For every a\\ , a reference non-polar structure is then ob- 
tained after a constrained relaxation within the PA/mbm 
group, until the maximum force and zz stress tensor com- 
ponents fall below 0.01 eV/A and 0.0001 eV/A^ respec- 
tively. Then, symmetry is broken displacing coherently 
the cations by hand, and a second relaxation is carried 



out without any imposed symmetry. 

For the (2/2) superlattice we have performed struc- 
tural relaxations under different in-plane strains, while 
keeping the global tetragonal symmetry. The misfit 
strain is defined as e = — — —, where an is our LDA 
theoretical lattice constant of cubic bulk STO (3.874 A). 

The dependence of the polarization (inferred from the 
bulk Born effective charges and the local atomic displace- 
ments), and the oxygen octahedra rotations and tiltings 
with the epitaxial strain can be seen in Fig. [2] The first 
conclusion that can be drawn is the existence of a strong 
coupling between FE distortions and AFD modes with 
strain. Both, the magnitude and the relevant directions 
(of the FE polarization and/or the rotation axis of oxy- 
gen octahedra) can be tuned by the amount of epitaxial 
strain. 

For large compressive strains the polarization along the 
z-direction (c-phase in Refs. [H, [laj), and the rotations 
of oxygen octahedra with respect the same z-axis are 
stabilized, while the tiltings are suppressed. Besides the 
effect of the epitaxial strain, the presence of the inter- 
faces plays a twofold role, as can be deduced comparing 
the superlattice and bulk results for pure PTO and STO 
in Fig. [5] (i) the magnitude of the polarization along 
z is homogeneous throughout the hetcrostructurc. The 
polarization mismatch at the interface is always smaller 
than 0.5 fiC/crn^, highlighting the large electrostatic cost 
of a polarization discontinuity between the layers 0, [l7j . 
The price to pay for poling the STO layer is a reduction 
in the polarization of PTO with respect to bulk [compare 
open and filled circles in Fig. [2ja)]. (ii) The magnitude 
of the angle of rotation is different for the different octa- 
hedra, a fact that points out a special coupling between 
the AFD modes and the FE polarization. It is remark- 
able that, although bulk PTO does not exhibit rotation 
of oxygen octahedra, the TiOg octahedra in the superlat- 
tice with PTO-likc environment inherits part of the AFD 
character of STO and arc forced to rotate. 

In the opposite limit, for large tensile strains, the polar- 
ization in the most stable configuration lies in the plane, 
along the [110] direction (aa-phase [l5, 16|). Note that, 
in this case, there is no electrostatic restriction to keep 
the in-plane polarization at the same value in the PTO 
and STO layers. Therefore, they arc decoupled and tend 
to the corresponding bulk values. The magnitude of the 
rotations along the z-axis are strongly reduced and also 
tend to bulk, while the tiltings along the [110] axis be- 
come the dominant AFD mode. As before, the tiltings 
in the PTO-likc octahedra are enhanced with respect the 
bulk values due to the interfacial coupling with the STO- 
like octahedra. 

Interestingly, at intermediate strains (around e ~ 0) 
the polarization rotates continuously from the c to the 
aa-phase (r-phase [HI, [IB])- Within this regime the elec- 
tromechanical response of the system {d^i and dn piezo- 
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electric constants) is enhanced, and at the STO lattice 
constant amounts to 0.30 nC/N, twice larger than the 
(^33 in Pb(Zro. 5X19.5)03 at room temperature [l^. The 
appearance of low symmetry phases where the polariza- 
tion is rotated away from the substrate normal has been 
already reported on PTO thin films grown on DySc03 
under tensile strains (+1.4%) [l^. However, we do ob- 
serve this rotation of polarization and the associated en- 
hanced piezoelectric response for much smaller strains. 
This is due to the increased stability of the r-phase in 
this system, which is caused by the constrain imposed 
by the STO on the polarization of the PTO layer. Our 
simulations suggest that, in order to reduce Pz, PTO 
prefers a rotation keeping large the magnitude of P over a 
monotonic reduction, a well known fact in FE perovskitc 

oxides fillii. 
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FIG. 2. (color online) (a) Polarization and (b) absolute value 
of the oxygen octahedra rotations and tiltings in monodomain 
(2/2) superlattices under different epitaxial strains. In (a), 
empty (patterned) circles correspond to bulk PTO (STO) un- 
der the same strain conditions. In (b), the rotation, (f), and 
tilting, 9, angles of the O octahedra (labeled as in Fig.[l]) are 
represented as triangles-up and triangles-down, respectively. 
As in (a), empty symbols refer to bulk values under strain. 
Top labels indicate strains induced by common substrates. 



As a summary of the coupling between strain and FE 
and AFD modes, we observe a rotation of the main axis 
of the system (defined by both the direction of the polar- 
ization and the rotation axis of the octahedra) from out- 



of-plane for compressive strains to in-plane for tensile- 
strains. For the AFD modes, this trend can be under- 
stood if we consider that the oxygen cage is very rigid. 
Then, as strains are applied, the system allows the TiOg 
octahedra to reorient to maintain their shape (see Fig. 3 
of Ref . ) . More importantly, here we observe an extra 
interface coupling between the FE polarization direction 
and AFD modes, that reveals itself in the different rota- 
tions of the octahedra types defined before. The largest 
difference is observed when strong compressive strains are 
applied: here the P"*" octahedra rotate more than PTO 
and STO ones, while P~ octahedra rotate much less. In 
order to understand this result let us discuss the origin of 
the coupling between FE and AFD modes in this system. 
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FIG. 3. (color online) (a) Diagram showing the change 
in distances between Pb and Sr ions with O anions at the 
PTO/STO interface under compressive strain for a P+ TiOe 
octahedra. (b) same as in (a) but for a P~ octahedra. Re- 
duction of distance and reinforcement of the Pb-O bond is 
shown by full red lines while an increase in the Pb-O distance 
and weakening of the bond is shown by dotted lines. Blue 
dotted lines represent Sr-O distances. Green arrows repre- 
sent out-of-plane displacement of A-cations (Pb or Sr), which 
magnitude is written in green. Yellow arrows represent the 
in-plane displacement of O. 



AFD distortions are usually regarded as purely steric 
phenomena, where the rotation of the octahedra takes 
place if the A-ion is small enough to let the B-O-B bond 
bend [2l|. A polar distortion of the A-cation along the 
positive z-axis reduces its distance with the oxygen ions 
of the Ti02 plane immediately above while increases the 
distance with the ones below [Fig. |31[a)]. Taking into 
account that the ionic radii for Sr^+ and Pb^+ are al- 
most the same (1.26 and 1.29 A, respectively) and using 
a steric model, we would expect that a stronger rotation 
would only be favored if the mean oxygen-A cation dis- 
tance increases, leaving space for the octahedra to rotate. 
However, we find that the mean A-0 distance is very sim- 
ilar for both P"'' and P~ octahedra, close to 2.78A , while 
the rotations are very different. In order to explain this 
difference in rotations, we propose that the driving force 
in the mixed AFD-FE-strain coupling in PTO/STO su- 
perlattices is of covalent nature, instead of steric. It is 
well known that a chemically active lone pair on the Pb 
ion, that allows for strong covalent hybridization with O, 
lies at the origin of FE in bulk PTO. For P+, due to 
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the coupling between FE and AFD distortions, not all 
the Pb-0 bonds are equivalent. In particular, having a 
shorter (2.45 A), and a longer one (2.93 A) allows the Pb 
complex to acquire a pseudo-tetrahedral shape typical of 
many covalent Pb^"*" compounds. For the P~ octahedra, 
the polar distortion increases the Pb-0 distances and the 
previous mechanism does not apply [Fig. [3jb)]. These 
results agree with recent ab-initio calculations that em- 
phasize the role of covalent interactions in the origin of 
AFD distortions 22] . As the in-plane strain is increased, 
the polarization rotates away from the z axis and this 
coupling is reduced making the in-plane rotations small 
and similar for all octahcdra-types when the values of the 
strain are larger than -1-1%. Under these tensile strains, 
the Pb displaces in-plane and both P+ and P~ become 
equivalent. 

We have also carried out simulations for different pe- 
riodicities, but fixing ay to the LDA theoretical one of 
STO. We perform both, unconstrained and constrained 
structural optimizations, where we impose a purely out- 
of-plane or in-plane polarization on the superlattice. Rel- 
ative energies and polarizations of the PTO and STO 
layers are gathered in Table U while the corresponding 
rotation angles can be found in the Supplemental mate- 
rials . The GS monodomain configuration displays 
both in-plane and out of plane polarizations, indepen- 
dently of n, although for n = 1 the GS r-phase reported 
in Table U is essentially degenerated with the c-phase (the 
difference in energy, 1.2 meV per 5 atom perovskite unit 
cell, is within the accuracy of our simulations). This del- 
icate competition was already observed by Bousquet et 
al, where the phonon frequency of the mode involving in- 
plane distortions in the (1/1) GS of Ref. was found to 
be of only 6 cm~^, close to become unstable. The small 
difference between the results in Table U and those in Ref. 
[lo| can be ascribed to small changes in the methodology. 
Larger periodicities of the superlattice seem to increase 
the range of stability of the r-region, as the difference in 
energies between this phase and the rest increases. For 
n > 2, within the PTO layer, P lies close to the diagonal 
of the perovskite unit cell (configuration labeled as [111] 
in Table HI . The GS can be considered as a condensa- 
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AFD2 + AFDj^y modes. In every 



case, Pz is nicely preserved at the PTO/STO interface, 
with a value between 30 and 35 /iC/cm^. in good agree- 
ment with previous first-principles simulations on clean 
PbTiOa/SrTiOs interfaces 

The structures obtained from the previous relaxations 
serve also as the starting point to answer the question 
about the evolution of the energy gain due to FE-AFD 
coupling with the periodicity of the superlattice. For this 
analysis we have focused on the [001] phase described 
above, since for this structure the separation of atomic 
displacements into polar distortions and oxygen octahe- 
dra rotations is trivial. Similar energy decompositions 
are expected for the GS [111] structure. The method to 
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Periodicity 

FIG. 4. (color online) Decomposition of the energy of the 
[001] phase into contributions from the polar modes, the AFD 
modes and the coupling between them. The reference energy 
corresponds to a structure where neither polar nor AFD in- 
stabilities are allowed to condense. 



disentangle the different energy contributions due to the 
polar displacements, oxygen octahedra rotations, and/or 
their coupling can be found in the Supplemental materi- 
als [23|. The results are depicted in Fig.|3]as a function of 
the superlattice periodicity. It immediately follows that 
for ultrashort periodicities the energy contribution from 
the coupling terms is negative and very large, increasing 
significantly the stability of the polarized structure. This 
points in the direction of previous theoretical results [lo| 
that support an "improper ferroelectricity" origin for the 
polarization in these ultrashort periodicities. However, 
the larger the periodicity the smaller the importance of 
the coupling term, that even changes its sign for n > 3. 
From this point on, we can consider our superlattices to 
behave as normal ferroelectrics, in good agreement with 
experimental results [according to Ref. H, (3/3) was the 
threshold periodicity above which the normal FE behav- 
ior was recovered]. 



In summary, our first-principles simulations show how 
the FE-AFD-strain coupling in monodomain PTO/STO 
superlattices produces a phase diagram much richer than 
initially envisaged. The driving force of the coupling is 
a combination of electrostatic and covalent effects. The 
new phases might contribute to the stabilization of the 
monodomain configuration over the recently observed 
and competing polydomain structures [2^. The experi- 
mental observation of the in-plane component of the po- 
larization in the superlattices remains to be confirmed. 
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CP-FP 228989-2 OxIDes. We acknowledge the computer 
resources, technical expertise and assistance provided by 
the Red Espanola de Supercomputacion. Calculations 
were also performed at the ATC group of the University 
of Cantabria. 
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TABLE I. Polarization and relative energies of tlie difTerent monodomain configurations for superlattices as a function of the 
periodicity of the supercell. In-plane strain corresponds to a STO substrate, with a theoretical in-plane lattice constant of oy 
= 3.874 A. GS stands for ground state. Polarizations in /^C/cm'^, energies in meV per five atom perovskite unit cell. 
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